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ABSTRACT 

Ellipsoidal red giant binaries are close binary systems where an unseen, relatively close companion distorts 
the red giant, leading to light variations as the red giant moves around its orbit. These binaries are likely to 
be the immediate evolutionary precursors of close binary planetary nebula and post-asymptotic giant branch 
and post-red giant branch stars. Due to the MACHO and OGLE photometric monitoring projects, the light 
variability nature of these ellipsoidal variables has been well studied. However, due to the lack of radial 
velocity curves, the nature of their masses, separations, and other orbital details has so far remained largely 
unknown. In order to improve this situation, we have carried out spectral monitoring observations of a large 
sample of 80 ellipsoidal variables in the Large Magellanic Cloud and we have derived radial velocity curves. At 
least 12 radial velocity points with good quality were obtained for most of the ellipsoidal variables. The radial 
velocity data are provided with this paper. Combining the photometric and radial velocity data, we present 
some statistical results related to the binary properties of these ellipsoidal variables. 

Subject headings: binaries: close - Magellanic Clouds - stars: AGB and post-AGB 


1. INTRODUCTION 

The variable red giants in the Large Magellanic Cloud 
(LMC) fall on six or more distinct sequences in a 
period-luminosity (PL) diagram, i.e., K-logP d i agram 
(iWood et all 119991: llta et all J2004t ISoszvnski et al.l 120071 : 
iFraser, Hawley, & Cookll20'08h . One of these sequences, se¬ 
quence E, consists of binary systems that are mainly red gian t 
ellipsoidal binaries (IWood et al.lll999t ISoszvhski et al.ll2004l) . 
In these ellipsoidal binary systems, the red giant is the pri¬ 
mary and it substantially fills its Roche lobe. The secondary, 
which is usually an unevolved main-sequence star, is unseen 
observationally in most cases. Due to the substantial filling of 
the Roche lobe, the red giant is distorted. Rotation of the dis¬ 
torted shape of the red giant as the binary progresses around 
its orbit causes a change in the apparent light seen by a distant 
observer. This leads to the characteristic light and velocity 
curves of ellipsoidal variables that have two cycles of light 
variation in one orbital period but only one cycle of radial ve¬ 
locity variation (e.g. INicholls et al1l2010h 
The sequence E stars are low-mass stars (m <1.85 M e ) or 
intermediate-mass stars (1.85 M 0 < m < 7.0 M Q ). They can 
lie on either the red giant branch (RGB) or the asymptotic gi¬ 
ant branch (AGB) and they are therefore in an evolutionary 
phase where the stellar radius is increasing. It is when the 
radius of the expanding red giant becomes a significant frac¬ 
tion of the binary separation that the ellipsoidal variability be¬ 
comes detectable with current surveys such as MACHO and 
OGLE. The orbital periods of sequence E stars in the LMC 
lie in the range ~30-1000 days and the light amplitude is usu¬ 
ally <0.3 mag in the MACHO red band Mr. Statistically, the 
sequence E stars make up approximately 0.5-2% of the RGB 
and AGB stars in the LMC. About 7% of the sequence E stars 
are eclipsing and about 10% of them have unusually shaped 
light curves which indicate significant eccentricity of the sys- 
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tern orbits (iSoszvhski et al.li2004l) . 

Due to the MACHO and OGLE projects, the variability 
nature of sequence E stars is well studied. However, at the 
present time, there are few radial velocity studies of the se¬ 
quence E stars, due _to the difficulties of long-time spectral 
monitoring. [Adams. Wood, & Cioni (2006) observed two se¬ 
quence E stars in the LMC, and their preliminary results in¬ 
dicated that the red giant is filling its Roche lobe and trans¬ 
ferring mass to the companion since the derive d mass o f the 
red giant is close to that of the red giant core. iNicholls et al.l 
(120101) carried out spectral obse r vatio ns on 11 s equence E 
stars (including the two in [Adams, Wood, & Cionil 12006 ) to 
derive radial velocity variations. The average full velocity am¬ 
plitude derived for those 11 sequence E samples was 43.4 km 
s' 1 , consistent with the inference that sequence E stars are red 
giants in binary systems with roughly solar-mass components. 

A problem with the current understanding of sequence E 
stars is that the eccentricity should be close to zero due to 
tidal interaction with the companion, y et about 10% of them 
have significant eccentricities (ISoszvhski et al.l 1200 4). Com¬ 
plete orbital solutions for ellipsoidal red giant binaries with a 
range of eccentricities could possibly show in which part of 
parameter spa ce eccentricity can be ma intained. As a start to 
such a study, INicholls & Woodl ( 2012 ) monitored the radial 
velocities of 7 eccentric sequence E stars and derived com¬ 
plete orbital solutions. 

In general, a knowledge of the complete set of orbital prop¬ 
erties (masses, separations, eccentricities, orientations) for a 
large sample of ellipsoidal red giant binaries in the LMC will 
provide a good resource for understanding the evolution of 
these close binary systems. These observed parameters can be 
used to constrain Monte Carlo simulations of the popu lation 
of sequence E stars, such as those in lNie, Wood, & Nichollsl 
(120120 . These Monte Carlo calculations provide estimates for 
the production rate of binary post-AGB stars, binary planetary 
nebulae (PNe), binary post-RGB stars, and luminous white 
dwarfs relative to the production rate of single-star post-AGB 
stars and PNe. The relative numbers of the various types of 
objects listed above, now quite well known in the LMC, can 
be used to parameterize the binary interaction process, espe- 
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daily common envelope evolution, and to calibrate models for 
the formation of cataclysmic variable stars, AM CVn systems, 
the progenitors of SNe la, etc. So knowledge of the parame¬ 
ters of a large sample of the precursor binaries should help us 
better understand the interaction processes and the evolution¬ 
ary fates of close red giant binaries. 

In this study, we present the results of radial velocity moni¬ 
toring of a large sample of ellipsoidal variables in the LMC. In 
Section[2] the design of our observing project and the data re¬ 
duction is described. In Section 0 we present observed prop¬ 
erties of binary stars based on the radial velocity and light 
variation data. Some statistical analyses are also provided and 
compared with theoretical models. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Selection of Objects 

We initially selected 8 6 sequ ence E candidates from those 
given in lSoszvriski et al.l C20041) . Only objects with I < 16.5 
mag were considered since good radial velocities could not be 
obtained in a reasonable time (20 minute exposure) for fainter 
objects. Also, the data quality of fainter stars is relatively poor 
in the MACHO and OGLE databases, so their light curves 
have large dispersions and this is not good for fitting orbital 
solution. 

All the sequence E candidates in ISoszynski et alJ (120041) 
and our 86 ellipsoidal variable candidates are plotted in the I— 
log P and VV'-log P planes in FigureQ] where W = / — 1.55 (V- 
I) is a reddening-free Weisenheit index ( Mador e 1982). It 
should also be noted that P is the orbital period, which is twice 
the period obtained by Fourier analysis of raw light curves 
since ellipsoidal variables have two maxima and two minima 
of the light curve per orbital period. The sequence E variables 
are often plotted in PL diagrams using the semi-period. 

Our aim was to select a sample of objects with as wide a 
set of parameters as possible. We therefore selected objects 
randomly throughout sequence E, as shown in FigureQ] sub¬ 
ject to the constraint I < 16.5. Note that the magnitude cut 
of 16.5 is about 2 mag below the RGB tip, so our sample in¬ 
cludes low-mass sequence E stars, even down to the masses of 
globular cluster red giants. We also made no attempt to select 
highly eccentric stars so the sampl e should include a range of 
eccentricities (note that lNicholls & Woodll2012l have already 
studied a sample of 7 eccentric sequence E stars). 

At a given orbital period, more massive binary systems will 
have larger separations so the red giant will need to get to a 
higher luminosity (radius) before the Roche lobe filling factor 
is large enough to produce detectable ellipsoidal light varia¬ 
tions. The more massive binaries should therefore preferen¬ 
tially lie in the high-luminosity side of se quence E in FigureQ] 
The results of INicholls & Woodl (120121) . shown in Figure Q] 
indicate that the higher-mass stars do indeed lie on the high- 
luminosity side of sequence E. In order to be sure of selecting 
a good sample of intermediate-mass stars (M >~1.85M Q ), 
we preferentially included extra stars above the blue line in 
FigureQ] Our final sample consists of 86 stars. 

With 86 ellipsoidal candidates, plus 11 sequen ce E stars 
whose radial velocities were obtained by INicholls et al.l 
(1201 Oh (green crosses in Figure [Tb and the 7 eccen tric se¬ 
quence E stars studied by INicholls & Woodl (1201 2l) (green 
squares in FigureQ] , we have a combined sample of approx¬ 
imately 100. The full sample will give us statistical informa¬ 
tion about properties such as masses, mass ratio, eccentricity, 
separation, Roche lobe filling factor, and red giant luminos- 



Fig. 1. — Location of our 86 ellipsoidal variable candidates (red dots) on 
the /-log P and VT-log P planes. The smal l black dots show the full sample 
of ellipsoidal variables in ISoszynski et alJ <20041) . The blue line in the top 
panel is a line parallel to the main part of sequence E. Stars above this line 
are brighter than normal for their period and are expected to be intermediate- 
mass stars. Objects in our sample that lie above the line are shown in both 
p lots as circled points. T he green squares are ellipsoidal variables studied by 
INicholls & Woodl <20121) and the area of the square is proportional to the mass 
of the red giant derived by them . The green crosses are ellipsoidal variables 
studied bvlNicholls et alJ <20101) . Due to the lack of V band data, 3 objects of 
INicholls et alJ (20101) are absent. 

ity for field red giants in the LMC. Derivation of some of 
these properties will require the use of a binary or bit m odeling 
tool such as the Wilson-Devinnev code (IWilson & Devinnevl 

fl97lt lWilsonl[l97lfl990l: IWilson et al.ll2009t). 

2.2. Spectral Observation 

The radial velocity observatio ns were taken using the 
Wide Field Spectrograph (WiFeS; iDopita et alJl2007l ifoloh 
mounted on the Australian National University 2.3m tele¬ 
scope at Siding Spring Observatory. WiFes is an integral 
field, double-beam, imaging-slicing spectrograph with a field 
of view 25x38 square arcsec, imaged onto 25 slits that are 1 
arcsec wide and 38 arcsec long. It has six gratings, giving 
high (R=7000) and low (R=3000) spectral resolutions. For 
our observations, the gratings B7000 (wavelength coverage 
of 4184-5580A) and 17000 (wavelength coverage of 6832- 
9120A) were chosen for the blue and red CCD, respectively. 
These two gratings give a two-pixel resolution R=7000, cor¬ 
responding to a 45 km s -1 velocity resolution. 

We carried out 18 weeks of radial velocity monitoring, from 
2010 September to 2012 March, to cover at least one orbital 
period for nearly all of the 86 objects. The observations were 
approximately evenly distributed throughout the 18 months, 
roughly one radial velocity observation per star per month. 
The exposure time was generally set to 300 s for objects with 
7~13 mag, 600 s for 7~14 mag, and 900 s for /—15 mag, 
which gives a signal-to-noise ratio (S/N) of at least 20. For 
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fainter objects (/ >16 mag), the exposure time was increased 
to 1200 s. We also increased the exposure time when observ¬ 
ing in bad weather, to guarantee a S/N of 20. 

For our observations, we chose the “stellar” mode exposure, 
in which case only 12 slits of the spectrograph were used. For 
flatfielding the QI-1 lamp was used and for wavelength cali¬ 
bration a Ne-Ar arc lamp exposure was taken at the beginning 
of the night. For velocity derivation, the radial velocity stan¬ 
dard star HR9014 was observed. In addition, the white dwarf 
star EG131 was observed so that telluric lines could be used 
to remove any zero point error in the wavelength calibration 
arising from spectrograph drift over the night. To achieve high 
S/N, the exposure times of HR9014 and EG 131 were set to 10 
s and 900 s, respectively. 

2.3. Data Reduction 
2.3.1. Spectrum Reduction 

The WiFes data reduction pipeline ( iDonita et al.ll2010l) was 
used for the spectrum reduction. The pipeline combines the 
calibration and science data and provides one-dimensional 
spectra with most of the cosmic rays and sky emission lines 
removed. Since all our objects are red, most of their flux is 
concentrated in the red region of the spectrum. Therefore, we 
did not reduce the blue beam spectra because these spectra 
are of low S/N. The red spectra (6832-9120A) contain many 
prominent telluric absorption lines as well as emission lines 
of water, OH, and O 2 . The pipeline-reduced spectra some¬ 
times contained some residual telluric lines and cosmic rays. 
Because of this, we checked by eye all the reduced spectra of 
the same object and, by comparison, removed residual cosmic 
rays and residual sky emission lines manually. 

Spectrally, red giant binaries with white dwarf or neu¬ 
tron s tar companions always show characteristic emission 
lines (fA~l lenlll984 t lKenvonlll986l : iMikolaiewska et aklll997t : 
iBelczvhski et al.ll2000t) . These objects are the symbiotic stars. 
In passing, we note that none of our ellipsoidal variables show 
emission lines in their spectra. This means that none of the 
companions are likely to be white dwarfs or neutron stars. 

2.3.2. Radial Velocity Calculation 

(1) . Relative radial velocity The relative radial velocity and 

its error were computed using the IRAF fxcor package. 
The fxcor package uses a Fourier cross-correlation 
method to find the wavelength shift between an object 
and a template spectrum in a specified cross-correlation 
region. The template was usually a radial velocity stan¬ 
dard star with well-determined radial velocity and a 
spectral type similar to that of the program object. In 
our case, all the 86 ellipsoidal variables are stars of K- 
M spectral type, so we choose as the template HR9014, 
which is a K5 star with a well-determined radial veloc¬ 
ity of -20.4 km s _1 . The cross correlation was made on 
the wavelength interval 8400-8750 A which contains 
the Ca II triplet lines and is relatively free of telluric 
lines. The heliocentric radial velocity (rzhetio(obj)), ob¬ 
tained from fxcor after heliocentric correction and the 
inclusion of the heliocentric radial velocity of HR9014, 
was saved along with its error. The error in Wheiio(obj) 
is normally below 4 km s _1 . If the error was larger than 
10 km s -1 (due to bad weather and low S/N), then the 
velocity point was excluded from our data set. 

(2) . Zero-point correction To check the velocity calibration, 

all spectra of each program star were cross correlated 


TABLE 1 

Radial velocities v t and 1<t Errors cr„ r in km s 1 . Stars are identified by 
their OGLE II R.A. and decl. 


OGLE 050659.79 -692540.4 

OGLU 051256.36 -684937.5 

HJD( 2450000+) 

v r 

°V, 

HJD(2450000+) 

V T 

OV, 

5463.06104 

246.14 

1.74 

5461.14209 

237.26 

1.70 

5513.99023 

251.84 

2.10 

5513.03125 

247.21 

1.64 

5563.00781 

284.43 

1.48 

5563.07324 

257.77 

1.54 

5590.98877 

270.53 

2.71 

5590.06592 

256.17 

2.27 

5603.96924 

257.47 

2.12 

5647.92578 

243.94 

1.51 

5649.00684 

232.40 

1.40 

5673.95508 

229.63 

2.68 

5678.97363 

264.23 

4.53 

5693.92676 

230.36 

1.53 

5679.01221 

250.28 

4.04 

5765.23242 

231.89 

2.03 

5764.26514 

246.01 

3.65 

5846.08301 

242.70 

1.92 

5783.19287 

239.24 

1.88 

5867.01807 

252.29 

1.39 

5848.15039 

273.73 

7.73 

5907.19336 

257.03 

1.75 

5868.03613 

282.79 

2.54 

5944.04883 

255.03 

1.35 

5906.11133 

270.63 

1.82 

5991.98584 

243.43 

1.50 

5944.17676 

242.22 

3.00 

- 

- 

- 


with a template consisting of a single spectrum of the 
telluric standard star EG131. EG131 is a white dwarf 
that radiates almost like a blackbody. Thus, its spec¬ 
trum has only telluric lines, making it an ideal template 
for the zero-point correction. In principle, the relative 
velocity of the telluric lines in different spectra should 
be zero, but, due to the movement of the CCD system 
and spectrograph between the taking of the arc spec¬ 
trum at the beginning of the night and the taking of 
the object spectrum, as well as the changes of the at¬ 
mospheric pressure, there can be shifts in velocity. To 
compute this zero-point velocity shift i/ zp (obj), program 
stars were cross correlated with EG131 in the wave¬ 
length interval of 8120-8370 A. This region is domi¬ 
nated by telluric lines and it is close to the region of 
Ca II triplet, so the zero-point correction in this re¬ 
gion should similar to that of Ca II triplet region. The 
zero-point velocity correction varied from about -10 to 
34 km s -1 across the many nights of the observation. 
HR9014 was also cross correlated with EG131 to ob¬ 
tain its zero point velocity correction (i/ zp (HR9014)). 

(3). Absolute radial velocity To compute the absolute value 
for the observed radial velocity v r including a zero- 
point correction, we use the formula 

v r = WhelioCobj) - V zp (obj) + H zp (HR9014) . (1) 

The final radial velocity data for 2 objects are given in Ta¬ 
ble!]] The full table is available online. 

3. RESULTS 

3.1. Properties of the Observed Ellipsoidal Variables 

Before presenting the properties of the observed ellipsoidal 
variables, we need to remove candidates that are not clearly 
ellipsoidal variables with the help of the observed radial ve¬ 
locity. Among our 86 objects, we found that 80 of them are 
real ellipsoidal variables, showing two light maxima and two 
minima but only one velocity maximum and minimum in one 
orbital period. The remaining 6 candidates do not clearly 
satisfy this requirement due to their low-velocity amplitude 
relative to the noise so they were removed. These 6 objects 
could be ellipsoidal variables whose orbital plane lies close 
to the plane of the sky or they could be sequence D stars 
as these stars lie close to the sequence E stars in the I -log P 
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and VV'-log P planes and they have small velocity amplitudes 
dNicholls et al.ll2009l) . For the remainder of this paper, we fo¬ 
cus on the 80 real ellipsoidal variables while we discuss the 6 
rejected objects in the Appendix. 

Table [2] shows various properties of the 80 ellipsoidal vari¬ 
ables. For all objects, the period given has_been re-derived 
because the original value from lSoszynski et aD (120041) is not 
accurate enough. The MACHO and OGLE II light curve data 
were taken more than 10 years before our radial velocity data, 
which was obtained between years 2010 and 2012. Because 
of the long time span between the two sets of data, a small 
error in the period can cause a significant error in the phase of 
the light curve projected forward by more than 10 years. To 
obtain a more reliable period, when OGLE III data existed, 
we combined OGLE III and OGLE II light curves (covering 
the interval 1998-2009) and used the phase dispersion mini¬ 
mization (PDM) method ( [Stellingwerfil 19781) to calculate the 
period from the combined light curve. If the OGLE III data 
was not available for an object, we combined the OGLE II 
and MACHO data, covering the interval 1992-2000. 

The effective temperature was calculated by co nverting 
(/ - K)o to r eff using spline fits to the data in lHoudashelt et all 
(I2000allbl) . The reason we used (/ - K) o as the color index 
is because most of our sequence E stars are K- or early- 
M-type stars, and for these red giants I - K varies much 
more with T e $ than J - K so that photometric errors are less 
important for I - K. The mean I magnitude was derived 
from the OGLE photometry, while the K magnitude was ob¬ 
tained from_the Two Micron All Sky Survey (2MASS) cat¬ 
alog dCutri et al. |2003|) . To remove reddening, we adopted 
E{B - V)=0.08 ( Keller & Woodll2006l) along with E(V — /) = 
1.38 x E(B - V) (ISchlegel, Finkbeiner, & Davisll 1 998l) . and 
E(V - K) = 2.744 x E(B - V) and A(K) = 0.35 x E(B - V) 
dRieke & Lebofskvill985l) . The bolometric correction BCjc 
was cal culated from (J-K)q using the data in lHoudashelt et all 
(I2000albl) and the bolometric luminosity was calculated from 
Ko and BCy, with a distant modulus of the LMC 18.54 
dKeller & Wood|[2006l) . 

Finally, in Table [2] we indicate whether a star is below (a 
“1” in the last column) or above (a “2” in the last column) the 
higher-mass line in Figure |T] The last column in Table [2] also 
pr ovides notes on unusua l variability characteristics as given 
in lSoszvriski et al.l (120041) . Stars for which there is OGLE III 
data available are also indicated. The data in Table[2]are avail¬ 
able online. 


3.2. Light Curves and Radial Velocity Curves 

We present samples of light and radial velocity curves for 
the ellipsoidal variables in Figure [2] Light and radial veloc¬ 
ity curves for all objects are available as online data. The 
time series data for the light curves are from the MACHO, 
OGLE II, and OGLE III databases if they are available. Each 
radial velocity curve has at least 12 good quality data points 
in one orbital period, good enough for a bi nary orbital solu¬ 
tion. As confirmed by iNicholls et al.l (|201 Of) , our observed 
ellipsoidal variables show two light maxima and two min¬ 
ima in one orbital period, while the velocity shows only one 
maximum and one minimum. Among our 80 objects, most 
of them (60/80=75%) have almost circular orbits, since they 
have equal light maxima and they have the same light curve 
widths for the first and second maxima (e.g., the top two ob¬ 
jects of Figure |2}. The other systems are eccentric binaries 
(20/80=25%) since they show different light maxima or light 
curve widths in one orbital period (e.g., the bottom two ob¬ 


jects of Figure [2). In addition, there are binary systems that 
contain a semi-regular pulsating red giant (8/80=10%) (e.g., 
o bject OGLE 052948.84 in Figure |2ji, similar to those noted 
by Nicholls & Woodl d2012l) . For these systems, both pulsa¬ 
tion theory and binary theory can be used to constrain and 
compare stellar parameters. 

3.3. The Velocity Amplitude 

The distribution of full radial velocity amplitude A RV for 
the 80 ellipso i dal var iable s plus the combined samp le from 
INicholls et al.l (120101) and INicholls & Woodl (120121) is pre¬ 
sented in the left panel of Figure 0 (red). From the figure, 
it can be seen that the velocity amplitude varies from ~5 to 
~ 130 km s -1 , with a peak near 30 km s -1 . Note that the num¬ 
ber of stars with velocity amplitudes larger than 80 km s _1 is 
small. 

Our observed velocity amplitude dis tributi on is also com- 
pared to the model prediction of iNie, Wood, & Nichollsl 
(120 12i) (the black dashed line in the left panel of Figure 3). 
The predicted distribution is for ellipsoidal red giant binaries 
on the top one magnitude of the RGB (870-2190 L 0 ) that have 
light amplitudes detectable by the OGLE II observations. Our 
comparison sub-sample of observations consists of 51 objects 
for which the luminosity lies on the top one magnitude of the 
RGB (the blue solid line in the left panel). Note that we com¬ 
pare to the model of INie, Wood. & N ichollsl d2012 ) whic h is 
calibrated on the OGLE II data of fSoszvhski et al.l ( 20041) be¬ 
cause that calibration is better than the one based on MACHO 
data due to the extra sensitivity of OGLE II observations to 
small-amplitude light variations. 

The cumulative distributions corresponding to the his¬ 
tograms of the ellipsoidal variables on the top one magnitude 
of the RGB in Figure [3] are shown in the right panel. A two- 
sample Kolmogorov-Smirnov test gives a probability of up to 
0.20 that the observation and model come from the same un¬ 
derlying distribution. There is thus a modest probability that 
the model is consistent with both the OGLE II photometry 
and the velocity amplitudes derived in this study. 


3.4. The Mass Function 

The binary mass function f\ of a binary system where the 
primary star (in our case, the red giant) is observable is 

K]P ml sin 3 i 

fy = T~F = t “ v> • ( 2 ) 

In G (nil + 

where K\ is the radial velocity semi-amplitude of the primary 
star, P is the orbital period, nt\ is the mass of the primary 
star, and m 2 is the mass of the secondary star. Note that f\ 
is a minimum estimate for m 2 . The distribution of f\ for our 
sample of red giant ellipsoidal variables in the LMC is shown 
in Figure H] The figure shows a peak at about 0.2 M Q , pre¬ 
sumably due to the dominant low-mass population, but there 
are a few objects that must be in more massive binaries where 
the secondary star is of intermediate- mass with m 2 > 4 M e . 
Red giants of thi s mass we re found in ellipsoidal binaries in 
the LMC bv INicholls & Woodl dMl . If we assume a ran¬ 
dom pole orientation for the binary orbit, which implies that 
the mean of sin 3 i = jk, and if we also assume that the mean 
mass ratio is q - 1, then for stars in the peak of the distribu¬ 
tion around f\ ~ 0.2 M G , the mean mass of the red giant is 
m\ ~ 1.35 M 0 . This is similar to peak in the mass distribution 
of red giants predicted by m odeling the star formation history 
of the LMC (e.g.. INie, Wood, & NicholIsilMl) . 
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TABLE 2 

Properties of the observed ellipsoidal variables 


No 

Object 

(OGLE II Name) 

P 

(day) 

/ 

(mag) 

V 

(mag) 

M b 

(mag) 

M r 

(mag) 

K 

(mag) 

A/ 

(mag) 

ARV 
(km s _I ) 

L 

(L 0 ) 

T’eff 

(K) 

R 

(Ro) 

Remark 

1 

050107.08-692036.9 

166.8 

14.05 

14.75 

- 

- 

13.14 

0.110 

70 

3278 

5773 

58 

2,ogle3 

2 

050222.40-691733.6 

107.8 

16.00 

17.53 

16.85 

15.87 

13.90 

0.070 

32 

581 

3990 

51 

1 

3 

050254.15-692013.8 

171.5 

15.67 

17.08 

- 

- 

13.82 

0.025 

35 

736 

4234 

51 

1 

4 

050258.71-684406.6 

433.8 

15.24 

16.96 

16.85 

15.65 

13.18 

0.040 

38 

1149 

4020 
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Note. — Columns 4—8 are mean magnitudes in different bands; Column 9 is the mean amplitude of light variability in I band; Column 10 is the full amplitude of radial velocity; 
Columns 11-13 are the luminosity, effective temperature and mean radius derived as described in the text; Column 14 contains remarks: number 1 for stars below the line in Figure ITI 
and number 2 for stars above the line, lower case letters stand for binary types: “e,” eccentric, “+sr” a binary whose red giant shows semi-regular variability as well as ellipsoidal 
variability, and “ogle3” denotes an object which has OGLE III light curve data. 
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Fig. 2.— Light and radial velocity curves of ellipsoidal variables. Mr and Mg are the MACHO red and blue magnitudes, respectively. The top two objects are 
ellipsoidal variables with circular orbits while the bottom two have eccentric orbits. For the I band data, chocolate points are from OGLE III and red points are 
from OGLE II. 
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Fig. 3.— Distribution of the full radial velocity amplitude ARV. Left panel: 
the red solid line shows the distribution for t he 8 0 ellipsoidal variables an d 
the combined sample of INicholls et al.l <20101) and INicholls & Wood <20121) , 
the blue solid line shows the observed distribution for 51 ellipsoidal vari¬ 
ables on the top one magnitude of the RGB, and the black dashed line shows 
the predicted distribution fo r ellipsoidal variables on the to p one magnitude 
of the RGB by the model of INie, Wood, & Nichollsl <20121) with the calibra¬ 
tion of OGLE II data. Right panel: the predicted cumulative distribution for 
ellipsoidal variables on the top one magnitude of the RGB. Bl ue is for obser¬ 
vation s from 51 objects, and black is for the model of Nie, Wood, & Nicholls 
<2012!) with the calibration of OGLE II data. Note that all the distribution 
peaks are normalized to unity. 



Fig. 4.— Binary mass function distribution for red giant ellipsoidal variables 
in the LMC. The red line den otes our obser v ation, the gr een line denotes 
the combined observations of Nicholls et al. (2010) and Nicholls & Wood 
(Ml, and the black dashed line is the sum of the two distributions. 


3.5. The PL Diagram 

The variations of light amplitude, velocity amplitude and 
mass function across the PL diagram are shown in Figure 0 
In the top panel, we can see that the higher-amplitude ellip¬ 
soidal variables tend to lie on the higher luminosity, shorter 
period _side of seguence E. This is just as demonstrated by 
ISoszynski et al.l (l2004l) and is presumably because as the red 
giant in a binary system (with a given orbital period) evolves 
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Fig. 5.— PL diagram for the our 80 ellipsoidal variables. The top panel 
shows the position on the PL diagram of objects with different full-amplitude 
A I, the middle panel shows the position of objects with different full- 
amplitude ARV and the bottom panel shows the position of objects with dif¬ 
ferent mass function fi. In all panels, filled circles den ote observations fro m 
this work and open triangles denote observations from Nicholls et al] 120101) 
and INic holls & W ood 120121) . 


to higher luminosity, it will expand to a greater filling fraction 
of its Roche lobe and hence be more distorted and have a large 
light variation amplitude. 

The middle panel of Figure [5] shows a tendency for higher- 
velocity amplitudes to be associated with shorter-period or¬ 
bits, as might be expected. Of course, low-velocity ampli¬ 
tudes can be seen for shorter- period orbits if the orbital plane 
is oriented close to the plane of the sky: there are a few objects 
that seem to fall in this category. 

The mass function also shows some correlation with posi¬ 
tion in the PL diagram (bottom panel of Figure[5]). The most 
massive objects, with f\ > 0.8 M e , lie on the upper part of 
sequence E in agreement with the results shown in FigureQ] 

3.6. Velocity Amplitude versus Light Amplitude 

In Figure [6j we present the light amplitude as a function 
of the velocity amplitude. Red symbols are our obser va¬ 
tion, a nd g reen symbols are observa tions from lNicholls et all 
(120lOll and INicholls & Woodl (120121) . The 20 eccentric ellip¬ 
soidal variables from our observatio ns and the 7 highly eccen¬ 
tric ellipsoidal variables studied by iNicholls & Woodl (120121) 
are marked as open triangles. As expected, there is generally 
no correlation of the light and velocity amplitudes. The for¬ 
mer depends on the Roche lobe filling factor (and inclination), 
while the latter depends on the stellar masses and the orbital 
separation. In principal, the orbital separation could influence 
the Roche lobe filling factor, but in practice the orbital separa¬ 
tion effectively determines the luminosity on the giant branch 
where the Roche lobe nearly fills and ellipsoidal variability 
becomes detectable. This suggests that the orbital separation 
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Fig. 6.— Relation between the velocity amplitude and light amplitude. Red: 
observations from this work; circles are for circula r binaries and t r iangle s are 
for eccentric bina r ies. G reen: observations from [Nicholls et al. <20100 and 
Nicholls & Wood (2012); circles are for circular binaries and triangles are 
for eccentric binaries. Note that some objects are overwritten because they 
have the same velocity and light curve amplitude. 
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Fig. 7.— Relation between velocity amplitude and luminosity. Symbols are 
as in Figure [6] 


(velocity amplitude) should depend on luminosity for red gi¬ 
ant ellipsoidal variables, but the light amplitude should not 
depend on orbital separation (velocity amplitude). 

3.7. Velocity Amplitude versus Luminosity 

As suggested in the last section, the velocity amplitude 
should show some correlation with luminosity. In Figure [7J 
we show the dependence of luminosity on velocity ampli¬ 
tude. The binaries with larger velocity amplitudes do tend to 
have lower luminosities although the effect is not prominent. 
This is as expected, since higher-velocity amplitude generally 
means a smaller separation and thus Roche lobe filling should 
occur at a lower luminosity. 

4. SUMMARY AND CONCLUSIONS 

We have presented radial velocity observations obtained for 
80 ellipsoidal red giant binaries in the LMC. This samp le 
is muc h la rger than the prev i ous sa mples of iNicholls et all 
(120101) and INicholls & Woodl (120121) . The mass function of 
the sample suggests that the typical mass of the red giants in 
the ellipsoidal binaries is ~1.35M 0 , in agreement with esti¬ 
mates derived from the star formation history of the LMC. 

The main purpose of this paper was to present radial veloc¬ 
ity data and basic observational properties of 80 ellipsoidal 
variables in the LMC. In future work, the radial velocity data 
will be combined with MACHO and OGLE photometric light 
curve data to give complete orbital solutions for these binary 
systems. This is possible because the distance, and hence 
the luminosity, of thes e LMC objects is well known (e.g. 
INicholls & Woodl 120121) . The results of the complete solu¬ 
tions will yield statistical distributions of masses, mass ra¬ 
tios, separations, and eccentricities in the period range ob¬ 
served. These distributions of orbital elements for binaries 
in the LMC can be compared with the sol ar vi cinit y statistic al 
data given in the classic pape r oflDuauennoy & Mavorl (119911) 
and the recent paper bv IRaghavan et al.l (120101) . It will be in¬ 
teresting to see if the same distributions exists in samples of 
binaries with different metallicity distributions and different 
star formation histories. 
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APPENDIX 

LIST OF OBJECTS REJECTED 

Here we present data for the 6 objects that we removed from the 86 ellipsoidal candidates due to their poor velocity data relative 
to the noise. Properties of these 6 stars are given in Table [3] and their light and velocity curves are presented in Figure [8] Their 
velocity data is given in Table [T] 
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TABLE 3 

Properties of Objects Rejected 
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0.200 

■ 

2298 

3517 

130 

l,ogle3 


OGLE 050716.07-69325 

9.7 P = 563.3 
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Fig. 8.— Light and radial velocity curves of the 6 objects in the Appendix. 
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